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We analyze systematic (classical) and fundamental (quan- 
tum) limitations of the sensitivity of optical magnetometers 
resulting from ac-Stark shifts. We show that in contrast to 
absorption-based techniques, the signal reduction associated 
with classical broadening can be compensated in magnetome- 
ters based on phase measurements using electromagnetically 
induced transparency (EIT). However due to ac-Stark asso- 
ciated quantum noise the signal-to-noise ratio of EIT-based 
magnetometers attains a maximum value at a certain laser in- 
tensity. This value is independent on the quantum statistics 
of the light and defines a standard quantum limit of sensitiv- 
ity. We demonstrate that an EIT-based optical magnetome- 
ter in Faraday configuration is the best candidate to achieve 
the highest sensitivity of magnetic field detection and give a 
detailed analysis of such a device. 



I. INTRODUCTION 

The detection of magnetic fields by optical means is a 
well developed technique with applications ranging from 
geology and medicine [Q to fundamental tests of viola- 
tions of parity and time-reversal symmetry 

In spite of their great variety, optical magnetometers 
can be divided in two basic classes. In the first class 
light absorption at a magnetic resonance is used to de- 
tect Zeeman level shifts, while the second class makes use 
of the associated changes of the index of refraction. So 
called optical pumping magnetometer (OPM) |lj as well 
as dark-state magnetometers based on absorption mea- 
surements belong to the first class. The recently de- 
veloped magnetometers based on phase-coherent atomic 
media [|]||] and the mean-field laser magnetometer of ref. 
H belong to the second class. 

If systematic measurement errors can be avoided, 
which in practice can be a challenging task, the smallest 
detectable Zeeman shift (in units of frequency) is deter- 
mined by the ratio of the noise level of the signal S to its 
rate of change with respect to frequency 



dS/du 



(1) 



A fundamental lower limit of S no i sc results from photon 
counting errors due to shot-noise of the probe electromag- 
netic wave. (dS/duj) , which characterizes a "quality 



factor" of the system, is determined by an effective width 
of the magnetic resonance. The ultimate goal of magne- 
tometer design is to minimize the noise level and the 
effective width at the same time. 

The width of magnetic resonances in optical magne- 
tometers is subject to two types of broadening: reso- 
nant power-broadening due to the coupling of the opti- 
cal fields to the probe-transition and a broadening due 
to ac-Stark shifts resulting from non-resonant couplings 
to other transitions. As shown in Q and || power- 
broadening limits the simultaneous minimization of noise 

and (d«Sydc<j) in absorption based magnetometers. In 
such devices increasing the probe laser power reduces the 
shot-noise but does reduce the signal at the same time. 
As a consequence the sensitivity saturates at a rather low 
power level. On the other hand, as shown in Q and ||, 
this effect can be compensated in a magnetometer that 
detects phase shifts of the probe electromagnetic wave 
propagating in an optically thick atomic medium under 
conditions of electromagnetically induced transparency 
(EIT) 0. Theoretically a complete elimination is possi- 
ble in a 3-level A-type system. 

In any real atomic system, however, there are non- 
resonant couplings to additional levels which lead to ac- 
Stark shifts and an additional broadening of the mag- 
netic resonance proportional to the laser intensity. In the 
present paper we analyze the influence of ac-Stark shifts 
and show that they (i) can diminish the magnetometer 
signal and (u) lead to additional noise contributions. We 
show that in absorption based devices ac-Stark broaden- 
ing leads to a further reduction of the signal. In contrast 
it only gives rise to a bias phase shift in an phase-sensitive 
EIT magnetometer. This bias phase shift can be cali- 
brated but is still a major source of systematic errors. It 
can be eliminated, if an EIT magnetometer with Faraday 
configuration is considered. 

However, in both, absorptive and dispersive type de- 
vices, ac-Stark shifts give also rise to fundamental noise 
contributions which increase with the laser power more 
rapidly than shot noise. Hence the magnetometer sen- 
sitivity decreases above a certain power level. The 
maximum value of sensitivity constitutes the standard 
quantum limit. For an EIT magnetometer based on 
phase-shift measurements this limit is determined by 
the dispersion-absorption ratio of the medium and the 
intensity-phase noise coupling due to the self-phase mod- 
ulation associated with ac-Stark shifts. 

We also discuss the possibility of further increasing the 
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sensitivity by means of non-classical light fields and show 
that the maximum sensitivity is essentially independent 
of the light statistics. 

The paper is organized as follows: In Sec. II we dis- 
cuss the fundamental broadening mechanisms of mag- 
netic resonances, power-broadening and ac-Stark associ- 
ated broadenings. It is shown in Sec. Ill that the classical 
signal reduction due to these broadenings can be compen- 
sated in phase-sensitive EIT magnetometers in contrast 
to absorption-based techniques. In Sec. IV fundamen- 
tal quantum noise sources are discussed and the stan- 
dard quantum limit of magnetometer sensitivity derived. 
A detailed analysis of an EIT-Faraday magnetometer is 
given in Sec. V and the prospects of using non-classical 
input states are discussed. 

II. BROADENING OF MAGNETIC 
RESONANCES 

Optical magnetometers measure in essence the position 
of certain resonances which are sensitive to magnetic level 
shifts. An important quantity that determines the signal 
strength of such a measurement is the width of the mag- 
netic resonance. As a rule the narrower the resonance, 
the easier it is to detect level shifts. 

Magnetic resonances with small natural width can be 
obtained e.g. by coupling Zeeman or hyperfine compo- 
nents of ground states in atoms either with an RF field 
or via an optical Raman transition. In an optical mag- 
netometer these ground-state sub-levels are then coupled 
by laser fields to excited atomic states. The optical cou- 
pling is also used to detect energy shifts of the ground- 
state sub-levels induced by a magnetic field. However, 
at the same time this coupling leads to a broadening of 
the magnetic resonances via two mechanisms: (i) power- 
broadening and (ii) broadening due to ac-Stark shifts. 

A. Power- broadening 

The first mechanism is power-broadening due to the 
resonant interaction with the probe transition. When 
the Rabi-frequency SI of the optical probe field exceeds 
the value 

«2t ~ VlTo, (2) 

where 70 is the unbroadened width of the magnetic res- 
onance and 7 the homogeneous linewidth of the opti- 
cal transition, the magnetic resonance becomes power- 
broadened. (Here and below we assume that 7 ^> 
7o-) The effective width scales linearly with the Rabi- 
frequency SI of the optical field or the square root of the 
corresponding power 

r efi = 70 + a + (3) 



a is some numerical pre-factor of order unity that de- 
pends on the specific model [§J9). This broadening effect 
leads to a substantial limitation of the signal in an optical 
pumping magnetometer, as shown in jH| and S. 

B. Broadening due to ac-Stark shifts 

The second broadening mechanism is due to non- 
resonant couplings of the probe electromagnetic wave 
with other than the probe transition and the associated 
ac-Stark shifts. The ac-Stark effect leads to a shift of the 
magnetic resonance of 

\si\ 2 

Aw ac _Stark = ~~ 7 — (4) 
^0 

where Ao is some effective detuning of non-resonant tran- 
sitions from the frequency of the probe field weighted 
with relative oscillator strengths. SI is again the Rabi- 
frequency of the probe field corresponding to the resonant 
probe transition. ( A is of course just a model-dependent 
coupling parameter. We have used this notation here for 
simplicity of the discussions.) 

In the classical limit and for a homogeneous laser in- 
tensity throughout the atomic vapor, there is only a con- 
stant frequency shift due to the ac-Stark effect. This shift 
can be calibrated. However, maximum signal is usually 
achieved when the atomic density is chosen such that 
there is a substantial absorption of the probe field. Hence 
when the probe Rabi-frequency exceeds the value 

n£> t ~ (5) 

the resonance frequency changes as a function of prop- 
agation through the medium. This leads to an effective 
inhomogeneous broadening of the magnetic resonance. 
For example, the transmission of a cell containing atoms 
with a Lorentzian magnetic resonance subject to ac-Stark 
shifts is determined by the integrated imaginary part of 
the susceptibility (x" = I m M) 

i dZX " (Z ^i dZ 7o 2 + (A + l^)IW (6) 

|Sl(z)| 2 characterizes the z-dependent power of the probe 
field and A the detuning from the un-shifted transition 
frequency. It is easy to see, that there is a broadening of 
the magnetic resonance depending on the magnitude of 
the ac-Stark shifts and the details of the absorption pro- 
cess. An important feature is that this broadening is pro- 
portional to the square of the Rabi-frequency or the laser 
power. Thus above a certain power level, determined by 
Eq. ([5]) ac-Stark associated broadening can exceed power 
broadening, which leads e.g. to further reduction of the 
signal in an optical pumping magnetometer. 
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III. COMPENSATION OF BROADENING 
EFFECTS IN EIT MAGNETOMETER 

We here demonstrate that the classical broadening 
mechanisms discussed in the previous section do not nec- 
essarily lead to a reduction of the magnetometer signal 
if phase measurement techniques are used. It has been 
shown in detail in [pi and J8| , that power- broadening can 
be completely compensated in a phase measurement by 
making use of EIT in optically dense A-type systems. 

The 3- level A configuration of an EIT magnetometer as 
well as the associated linear susceptibility spectrum of the 
probe field are shown in Fig.l. Here and in the following 
we consider closed systems i.e. we assume that there 
are no effective decay mechanism due to time-of-flight 
limitations. The upper level of the probe-field transition 
| a) «-> | b) is coupled to a meta-stable lower level |c) by a 
coherent and strong driving field of Rabi- frequency fi<j- 
The probe field Rabi-frequency is denoted as O p (f2 p -C 
fid) and the coherence decay rate of the probe transition 
as 7. A is the one-photon detuning of the drive field 
and S the two-photon detuning. The transverse decay 
rate of the two-photon resonance (magnetic resonance) 
is denoted as 70. It is assumed that the corresponding 
population exchange between the ground-state sub-levels 
is small and will be neglected 

As in the case of an OPM there is power-broadening 
in an EIT magnetometer as soon as > ^770. A 
unique property of an EIT resonance is however that 
the dispersion-absorption ratio of the optical transition 
is given by the inverse of the width of the ground-state 
transition 70 and is independent on the drive power if 
\&d\ > v / 77o- Under conditions of one-photon resonance 
(A = 0) one finds for small two-photon detuning 

x ' = Re[ x ]~--— *- , (7) 

Pd\ 2 +77o 

X" = Im \x] ~ m | 2 ° • ( 8 ) 

The residual absorption at the EIT resonance decreases 
with increasing laser power in the same way as the dis- 
persion. Thus in a phase shift measurement power broad- 
ening can be compensated by increasing the density and 
keeping a constant optical depth of the medium. 

Similarly one finds that as long as the drive-field Rabi- 
frequency is large compared to probe-induced ac-Stark 
shifts, which is very well satisfied, ac-Stark shifts of the 
magnetic resonance (eq.(^)) lead only to a bias phase 
shift. 

^ac-Stark ^ Az^f^, (9) 

Jo A o 

where L is the length of the atomic vapor cell. This 
phase shift can in principle be calibrated but gives rise 
to systematic errors. As will be discussed in detail later 
on, there is no such bias phase shift in a resonant Faraday 
configuration. 



We conclude this section by emphasizing that in phase- 
detection schemes based on EIT the detrimental (classi- 
cal) effects of power-broadening and ac-Stark associated 
broadening are eliminated. In the following section we 
will discuss the fundamental quantum noise sources of 
such magnetometer schemes. 

IV. QUANTUM-NOISE LIMIT OF MAGNETIC 
FIELD DETECTION VIA OPTICAL PHASE 
SHIFTS IN THE PRESENCE OF AC-STARK 
EFFECTS 

The problem of sensitive detection of phase shifts is 
common in optics. On the quantum level, the sensitivity 
of such kind of measurements is restricted by (i) vacuum 
fluctuations in the system and (ii) self-action noise due 
to nonlinearities in the system, as for example caused by 
ac-Stark shifts. The simultaneous presence of both noises 
usually leads to an absolute limit of the sensitivity. 

Let us discuss this problem for the particular case 
of optical magnetometry based on phase-shift measure- 
ments in an atomic medium. The ultimate limit for the 
smallest detectable phase shift is set by the generalized 
uncertainty relation between phase- A<fi = <j> — (4>) 
[ p"l| and photon- number fluctuations An = n — (n) of the 
output field. 

(A^ 2 )(An 2 )>l + i({A0,An}) 2 , (10) 

where {,} denotes the anti-commutator. If phase- and 
photon-number fluctuations are uncorrelated, the second 
term on the r.h.s. vanishes and one recovers the familiar 
Heisenberg relation. In any real magnetometer schemes 
phase and intensity fluctuations are however correlated 
due to e.g. ac-Stark shifts (self phase modulation), and 
thus the second term in Eq. (|lOj) is in general nonzero. 
When the intensity-phase coupling is small, it can be 
characterized by a linear coupling coefficient (3 in the 
form A</> = A(/>o+/3An, where A(/>o denotes phase fluctu- 
ations not correlated to intensity fluctuations. Thus we 
find 

<A0 2 )>^y+/3W). (11) 

The signal phase accumulated during the propagation 
through an atomic vapor cell is proportional to the Zee- 
man splitting Aujb, the length of the cell L, and the 
dispersion of the real part of the susceptibility at the 
laser frequency dx'/duj. The cell length is restricted by 
the absorption at the laser frequency, and a reasonable 
upper limit for L is the (amplitude) absorption length 

L ahs = (7r X "/\y 1 . 

Thus the maximum phase shift is 
1 dv' 

A0| max = — -f- Alu b . (12) 
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One recognizes, that the sensitivity of phase measure- 
ments to Zeeman shifts is determined by the dispersion- 
absorption ratio (l/x") dx'/duj. 

The limit for the smallest detectable Zeeman shift is 
therefore given by 



Auj b \ 





-l 







(An 2 



p 2 (An 2 



1/2 



(13) 



Under the condition, that the dispersion-absorption ra- 
tio is independent on the laser power, the r.h.s. of this 
expression is minimized when (An 2 )| opt = /3 . There- 
fore there is an absolute lower limit or "quantum limit" 
of magnetic field detection via phase-shift measurements 
independent on the photon-number fluctuations 



Aw, 



X 



i_pV 



2/3. 



(14) 



The absorption-dispersion ratio of a magnetic resonance 
is usually given by its natural width, which can be rather 
small if a two-photon Raman process between Zeeman- 
or hyperfine components is used as in an EIT magne- 
tometer. 

We will show later on that different measurement 
strategies as well as the use of non-classical light fields 
do in general not improve this result. 



V. EIT-BASED FARADAY MAGNETOMETER 

Let us now discuss in detail an EIT magnetometer in 
resonant nonlinear Faraday configuration. For this we 
consider the propagation of a strong, linear polarized 
light field through an optically dense medium, consist- 
ing of resonant A-type systems (atoms, quantum wells 
etc.) as shown in Fig. 2. For simplicity we ignore opti- 
cal pumping into lower states other than those shown in 
the figure and assume a closed system. For a resonant 
J = 1 — > J = transition (say), optical pumping into 
the lower mj = state depletes both states mj = ±1 
in the same way and thus effectively diminishes the op- 
tical density but does not affect the signal. Symmetric 
re-pumping can be used to maintain the population in 
the relevant sub-system without affecting the detection 
scheme. We include a dephasing of the ground-state co- 
herence with rate 70 and a population exchange rate be- 
tween the ground states 7o r . 

The two circular components E— and E+ of the linear 
polarized light generate a coherent superposition (dark 
state) of the states \b±) = \ J = l,mj = ±1). A mag- 
netic field parallel to the propagation axis leads to an 
anti-symmetric level shift of \b±) and thus by virtue of 
the large linear dispersion at an EIT-resonance to an op- 
posite change in the index of refraction for both compo- 
nents. As a result the polarization direction is rotated, 
which is the so-called resonant nonlinear Faraday effect 
Jl2[. The difference to the linear Faraday effect is the 



presence of the intensity-dependent dark resonance gen- 
erated by the action of the strong laser field as opposed to 
a usual absorption resonance in the weak-field limit. The 
rotation of the plane of polarization at the output can be 
measured by detecting the intensity difference of two lin- 
ear polarized components ±45° rotated with respect to 
the input polarization. 

An aspect of the system, which becomes particularly 
important when strong fields are considered, are non- 
resonant couplings of the two circular components to 
other levels, which to lowest order give rise to ac-Stark 
shifts of the states \b±). In a Faraday configuration the 
ac-Stark shifts of |& + ) and |&_) are exactly equal and 
opposite in sign due to symmetry and thus there is no 
average effect on the signal and no bias phase shift or 
rotation. Thus the Faraday magnetometer is not sub- 
ject to systematic errors associated with ac-Stark shifts. 
However, as mentioned before, ac-Stark shifts cause a 
coupling between intensity and phase fluctuations which 
need to be taken into account. 



A. Detection scheme 

We here consider the detection scheme shown in Fig. 3. 
A strong linear polarized field initially polarized in x di- 
rection propagates through a cell of length L with the 
magneto-optic medium. Due to the nonlinear Faraday 
effect the plane of polarization is rotated by an angle 
0/2. 

In order to detect this angle the intensity difference of 
the two orthogonal output directions 1 and 2 is measured. 
The operator for the number of counts is given by 

n = C^dt (i£{t)E+(t) - E^{t)E+(t)). (15) 

where E^ denote the positive and negative frequency 
part of the output electric field operators, t m is 
the measurement time, and C — 2eocA/%VQ, A be- 
ing the beam cross-section and vq the resonance fre- 
quency. Making use of the field commutation relations 
[E+ 2 (L,t),E^ 2 (L,t')} = C- x 5(t-t') and = 0, 

we can express the mean number of counts as well as the 
fluctuations in terms of normal-ordered correlation func- 
tions. The latter allows to apply a c-number approach 
where the operators E are approximated by stochastic 
complex functions E. 



{n) = {n 2 ) - {nx) 

-2\ _ /A„2\ 



(An 2 ) = (An 2 ) + (m) + (n 2 ). 



(16) 
(17) 



where 711,2 follows form Eq.(|15|) by replacing the field 
operators by c-numbers 



™1,2 



= C 



J^dtE^ 2 (L,t)E+ 2 (L,t). 



(18) 
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In the usual configuration only the x-polarized com- 
ponent of the input field is excited and we will restrict 
the discussion to a vacuum input of the y-polarized 
component. The propagation of the field through the 
magneto-optical medium is most conveniently described 
in terms of right and left circular components E± — 
(1/^/2){E x ± iEy^J, and we therefore have 



obtained from the stationary solution of the c-number 
Bloch equations for the atomic populations 



O&- b- 



06+ 6+ 



= -iC J dt^EZE^ - E+E±y (19) and polarizations 



— 70r(oh-6- - CT6+6+) + 7i' cr aa 
-i(fi_(7 Q 6- - C.C.), 
70r(cb- 6- - CT6+6+) + 7rO- Q a 
-i(ft + (T ab+ - C.C.), 



The propagation of the circular components can be char- 
acterized by two parameters, the intensity transmission 
coefficient rj and the phase shift (j)±(L, t) of the respective 
component at the output. 



E+(L,t)=E+(0,t)^e^ L ' t \ 



(20) 



In the limit of small magnetic fields the absorption of 
both circular components is identical for symmetry rea- 
sons i.e. there is no dichroism. With this we obtain for 
cw-input fields 



(n) = r)(n x } ia sin(f> sig ps rj(n x ) in (j>. 



(21) 



where s ; g = </>+(£) — (j)-(L) is the (stationary) signal 
phase shift. Similarly we can estimate the fluctuations in 
lowest order of the small rotation angle 4> in the case of 
an initially coherent field 



(An ) = n(n x )in + r\ (n x )( 



(22) 



The first term corresponds to the vacuum noise level and 
the second term proportional to 

w 2 > = // dtdt ' (m),m')) (23) 

describes fluctuations due to an intensity-phase noise 
coupling in the medium, ((a, b) = ([a — {a)][b — (&)])) 

In the following we calculate the loss factor r], the sig- 
nal phase shift </> s ; g and the fluctuations (6(j) 2 } due to the 
intensity-phase noise coupling for the EIT-Faraday mag- 
netometer. 



B. Medium susceptibility and field propagation 

The stationary propagation of the right and left cir- 
cular polarized electric field components through the 
atomic vapor is described by Maxwell equations in 
slowly-varying amplitude and phase approximation 



d , ii>Q 



(24) 



N is the atomic number density, p± are the dipole mo- 
ments of the respective transitions, and <Jb± a are the 
c-number analogues of the atomic lowering operators 
o&±a = |6±)(a|. Analytic expressions for CT6±a can be 



cr a b± 



0&- 6+ 

where 



~ra6±0- o 6± — iQ±{(T b ±b± — Oaa) 
— if2*p06+ 6±, 

-r b -6+0-&-6+ - iQ-<7 a b + + i^+O-6-a, 



r a &± = 7 - 



— +l ^A + S ± ±- 

IV 6+ = 7o + 7or + i (So + S+ - S-) 



(25) 
(26) 



(27) 
(28) 



7 r is the radiative linewidth of the transitions |o) — > \b±), 
and 7 is the homogeneous transverse linewidth of the op- 
tical transitions \a) — > \b±). So is the Zeeman splitting 
and 5± are the ac-Stark shifts of levels \b±). Cl± are 
the complex Rabi-frequencies of the two optical fields, 
il± = p±E±/h. We have disregarded Langevin noise 
forces in Eqs.([25[-p6|) associated with spontaneous emis- 
sion and collisional decay processes, since it was shown 
in H that atomic noises have a negligible effect on the 
magnetometer sensitivity. 

We calculate the stationary solutions of the Bloch- 
equations by considering only the lowest order in 70, 7o r , 
Sq and S±. In this limit we find 



o- a 6 ± 



«f2±(7o|fi 



+ 



70] 



l^±| 2 ) 



-|n| a ) 
2n±\n 



|fi| 2 (27(27 0r +7o) 
5o\ 

2 J|fi| 2 (2 7 (27 0r + 70 ) + |r!| 2 ) 

a n± hord^+l 4 + 1«_| 4 ) + 2 7o |r! T | 2 |r!| 2 



+ 



(29) 



7 



|r>| 4 (27(2 7 o r + 7o) + |0| 2 ) 



where \fl(z)\ 2 = \fl-(z)\ 2 + \n+(z)\ 2 . Usually the co- 
herence decay between the ground levels dominates the 
population exchange and thus 70 3> 7o r - 

It is convenient to separately consider the spatial evo- 
lution of amplitudes and phases of the complex Rabi- 
frequencies £l±(z) = \Q±(z)\e z ^ ± ^ . The intensities of 
the two fields are attenuated in the same way 



707r 



\n+\ 2 \n\ 



\n\ 2 (2 7 o7 + W 2 )' 



(30) 



where k = (3/4ir)N\ 2 . 

Eq. ( |30| ) can be easily solved when the length L of the 
cell is small enough, such that |f2(L)| 2 ^> 2770. In the 
Faraday set-up discussed here S1±(0) = Q(0)/V2, and 
therefore |f2±(z)| 2 = |f2(z)| 2 /2. We thus arrive at 
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|tt(z)| 2 = |tt(0)| 2 ( 



707 r KZ 



2\m\ 2 

|r!(0)| 2 (l-a z 



(31) 



It is interesting to note that under conditions of EIT the 
residual absorption is not exponential but linear. The 
intensity transmission coefficient is then given by 



V 



1 — UqL 



(32) 



The approximation |fi(L)| 2 3> 2770 sets an upper limit 
for the losses, such that 1 > 77 ~» 277 /|fi(0)| 2 . 
Similarly we find the phase equations 



d_ = K 7r A 7o + 7(<W2-(5_) 

dz ~ 2 7 2 70 7+|fi| 2 

d_ = K 7 r A70- l(5 /2 + S + ) 

dz + 2 7 2 70 7+|fi| 2 



(33) 
(34) 



The contributions from the one-photon detuning A can- 
cel when the relative phase <fi — 0+ — 4>- is considered 



d 

Tz 1 



KJ 

2 



n\ 2 



Ifil 2 



(35) 



The first term describes the signal-phase shift due to a 
magnetic field and the second term the ac-Stark contri- 
bution. Integration of Eq. (p5|) yields for the signal 



So 

7o 



In 



n(o) 



n(L) 



and the ac-Stark contribution 



«7r 



dz 



5+(z,t) - 5-(z,t) 

rap ■ 



(36) 



(37) 



C. Ac-Stark shifts and associated noise 

Let us now discuss the average ac-Stark shift and the 
corresponding noise contributions. For this we first con- 
sider the effect of an off-resonant quantized field on the 
energy of a single atom in lowest non vanishing order of 
perturbation. We then generalize the results for the aver- 
age ac-Stark shift and its fluctuations to an ensemble of 
atoms by making the physically reasonable assumption 
that ac-Stark shifts of different atoms are uncorrelated. 

We decompose the Hamiltonian of the single atom in- 
teracting with the quantized field in a rotating frame in 
the form H = Hq + H$, where Hq is the unperturbed 
part 



H Q = H™ d + hLo b _ b+ \b_)(b_\ 

+hA ab+ \a)(a\+h^2A j \c :j 



(38) 



A a b + = Lu ab+ — vq and Aj — uJ Cj b + — v§ are the detunings 
of the I a) — \b+) and \cj) — \b+) transitions. 

H s = -p-\a){b-\E±-p + \a){b+\E+ 

-^(pi + \<?){K\H + Pi-\<h)(t>-\*£) + ad J- (39) 

j 

describes the resonant and non-resonant couplings of the 
quantized fields to the atom. The non-resonant couplings 
to the excited states \cj) cause ac-Stark shifts. We here 
have assumed that both fields are nearly monochromatic 
and have set the energy of level \b+) equal to zero. p 3 ■ ± 



are the dipole moments of the transitions 



We proceed by formally eliminating the excited states 
\cj) by means of a canonical transformation in second 
order perturbation 

H = exp(S) H cxp(-S*) ~ H + [S, H] + [S, [S, H}} , (40) 
where S obeys the equation 

[S, H ] = {Pi + \ c i) ( b +\$+ + Pi-\ci) (b- \E± + adj.) 

j 

(41) 

Under conditions of exact two-photon resonance for the 
fields we obtain the transformation operator 



S = J2 



pj 



A |Cj >{6 + |^+ + ^_| Cj )(6_| J B± _ adj 



(42) 



Assuming that the population of all excited levels is 
small, we eventually find for the transformed Hamilto- 
nian 

H ~ H - p + \a)(b + \E+ - p_\a)(b_\E± 



b-)(b-\EIE^_ 



E 



hA 3 



b+)(b-\E+E± + \b-)(b+\EZE^ 



(43) 



Let us assume now, that Aj is much larger than the 
natural width of the excited states, and therefore the 
population transfer due to the non-resonant coupling is 
negligible. We identify J2j p]±/Aj — ► p 2 /Ao, where Ao 
is some effective detuning. The dipole moments pj+ and 
pj- have usually alternating signs for different excited 
states |cj). We therefore set Y^j Pj+Pj- /Ai = 0. Then 
the ac-Stark shift of the single atom can be represented 
by the operator expression 



Si 



(t) 



P" 



h 2 A 



E 



± (z u i)E^(z h i), 



(44) 



where I specifies the atom and zi its location. Thus we 
find for the average ac-Stark shift 



G 



(5 l ± (t)) = -^-(E ± (z h t)E+( Zll t)) 



M*i,t)\' 

2An 



, (45) 



where p 2 \(E ± {z u t))\ 2 /h 2 = p 2 \(E(z h t))\ 2 /2h 2 = 
\Sl(zi, t)\ 2 /2. Similarly we obtain for the second- 
order moments of the ac-Stark shifts (x, y) = 

((x-(x))(y-(y))} 



(8 l ± (t)6 l ± (t')) 

4 

= ^{Egizt^E+izt^Egizut^E+izut')), (46) 



(S l + (t)S l _(t')) 

= ^(E+izu^E+iz^EZizut'^izut')), (47) 

or after normal ordering 
(5 l ± (t)5 l ± (t')) = 

[(Eg ( Zl ,t)E+{z u t)E ± (z h t')E±(z h t')) 
6{t - t') 



h 4 A 2 



C 

(6 l + (t)S l _(t>)) = 



-{E±{z u t)EZ{ Zl ,t)) 



(48) 



(E7 (zi , t)E^ (zi , t)EZ (zi , t')E±(zi ,t')) . (49) 



h 4 A 2 



The first terms in Eqs.(|48|) and (|49j) correspond to classi- 
cal fluctuations, while the second term in (jig) is vacuum 
or shot noise. If the applied fields are in a coherent state 
only the shot noise term survives. In any practical re- 
alizations there are however large excess noise contribu- 
tions and the first terms are usually the dominant ones. 
We will show that all excess noise contributions are can- 
celed in a Faraday magnetometer and only the vacuum 
contribution survives. 

We generalize the above single-atom results to an en- 
semble of atoms assuming independent fluctuations of the 
ac-Stark shifts of different atoms, i.e. 



Ojk, 



(50) 



where {/i, v} 6 {+,—}. We introduce the continuous 
variable 



5 ± (z,t)=Lj2^-^i(t)- 



(51) 



(S±(z,t),S ± (z',t')) 



YaJ 6{z - z) 



(E±(z,t)E+{z,t),E±{z,t!)E+{z,t')) 
S(t - 1') 



C 



-(EZ(z,t)E+(z,t)) 



(53) 



and 



(MM),M^')> 



Lp A 
h*A 2 



5{z-z') 



(E±(z, t)E+(z, t), E±(z, t')E+(z, t')) 



(54) 



We here have used that in continuum approximation for 
any smooth function f(z) holds 

LY,S(z- Zl )5(z' - zi)f(zt) = S(z - z')f(z). (55) 



It is now straight forward to evaluate the quadratic 
deviation of the relative ac-Stark shift 

/ / 5 + (z, t) - L(z, t)\ ( h{z', t') - L (z', t') \ \ 



m{z)\ 2 

S(z- z')5(t-t') 



2\n(z'W 



p 2 L 



2h 2 CA 2 \n(z)\ 2 



(56) 



We note that the classical excess noise contributions ex- 
actly cancel and only the vacuum contribution is left over. 
Due to the intrinsic balancing in the EIT-Faraday mag- 
netometer excess noise contributions are automatically 
canceled. This is an important advantage of the Fara- 
day configuration as compared to the asymmetric EIT- 
magnetometer discussed in fl a nd ||. 

Using Eqs.(^3|), ( |37| ) and (|56|) we eventually find for 
the phase fluctuations due to ac-Stark shifts 



1 



1 ^ltp 2 L 
tm 4Ag h 2 C Jo \tt(z 



-dz 



(57) 



D. Signal-to- noise ratio and minimum detectable 
Zeeman shift 



The minimum detectable Zeeman shift is obtained by 
setting the mean number of counts 



In a continuum approximation, ^ - — > (1/L) J L dz, and 
we have 



<4(M)> = t?- 

Similarly 



(EZ{z,t)Et{z,t)) = 



2A 



(52) 



(h) = T] (n x ) in <p sig = -T) (n x ) in —In [rj *] (58) 

7o 

equal to the quantum mechanical uncertainty 

1/2 



(A: 



■2,1/2 = 



r/(n x ) in + r/ 2 (n x ) 2 n (6(/> 2 } 

, \m\ 4 



■•q (1 — rj) ln(?y 1 ) 



1/2 



(59) 
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This yields the signal-to-noise ratio 



SNR = 



S 2 2 - n 1/2 

-4(n x )i n ?j In 2 (77 x ) 

7o 

\m\ 4 

Ag 7o 2 



(60) 



which is maximized for an optimal power of the field cor- 
responding to 



M0)\l pt 



T) (1 - 7?) 01(77 r ) 



A 7o- 



(61) 



Substituting the optimum Rabi-frequency ( |6l| ) into j60| ) 
yields a maximum SNR for 77 « 0.06. Thus we find the 
quantum limit for the detection of Zeeman level shifts 



where 



<5 SQL =7o/ 



/ 



7r 



3 A 2 1 



A 8tt A j t r , 



1/2 



77 ln^-i) 



1/4 



(62) 



(63) 



is a numerical factor which varies between 1 and 2 for 
r\ = 0.01 ■ • • 0.8. (Note that 77 is the transmission coeffi- 
cient under conditions of EIT. Without EIT the medium 
would be totally opaque.) In Fig. 4 we have shown the 
minimum detectable Zeeman splitting (proportional to 
the magnetic field) as function of the laser input power 
for different transmission coefficients. 

One clearly sees that for small laser powers shot-noise 
is dominant, while for larger laser powers ac-Stark associ- 
ated fluctuations take over. Also shown is the saturation 
behavior of an OPM g. Due to power broadening the 
sensitivity of an OPM saturates as soon as the Rabi- 
frequency reaches the value y / 77o. In the EIT-Faraday 
magnetometer, on the other hand, the optimum Rabi- 
frequency corresponding to the quantum limit is of the 
order of y'Ao7o- Since Ao S> 7 much higher sensitivities 
can be achieved here. 



E. Compensation of ac-Stark associated noise by use 
of non-classical input fields 



In the case of an optical magnetometer, ac-Stark shifts 
appear due to non-resonant nonlinearities and it would 
seem that these shifts can in principle be compensated 
by an adapted measurement strategy and the use of non- 
classical light. An essential condition for such methods is 
however that the system is nearly lossless in order to pre- 
serve the non-classical state of light. On the other hand, 
as discussed above, the maximum signal in an optical 
magnetometer is achieved under conditions of substan- 
tial absorption. (We note that the SNR is proportional 
to ln(rj) 2 .) We will show in the following with simple es- 
timates that this feature makes it impossible to increase 
the sensitivity by using non-classical light. 

Let us consider the simplest example of compensation 
of ac-Stark associated noise by non-classical light. We 
assume, that the slowly varying field operators in the 
Heisenberg picture are represented in the form E± = 
(E) + e±, where e± is the fluctuation part. To discuss 
the compensation of ac-Stark effects let us disregard the 
resonant coupling with the medium and the associated 
absorption. Then we find that the field fluctuations at 
the end of the vapor cell can be written as 



e-{t,L) 
e+{t,L) 



e_(f,0)+i 
= e+(t, 0)-i 



K"f T L 



2A 

K~f T L " 



e_(t,0) + e±(f,0) 



2A 



e+(t,0)+e+(t,0) 



(64) 
(65) 



The second terms in these equations are due to ac-Stark 
shifts. One can see that the uncertainty of the phase 
difference increases as a result of ac-Stark shifts, which 
leads to the sensitivity restriction, discussed above. 

Let us assume now, that the incident field is squeezed 
in such a way, that the operators of the field fluctuations 
at the input obey the relations 



e_(t,0) = 
e+(t,0) 



.kj t L 
7tt^ — (e 
2A v 



et) 



2A 



(e+ + e+). 



(66) 
(67) 



Here e± are free-field operators (the corresponding state 
is the field vacuum), which obey the commutation rela- 
tions [ei(7j),e±(i')] = C-^it-t') and [g|(t),g|(f )] = 0. 
Then, in the absence of losses, the effects of ac-Stark 
shifts are completely compensated in the output and the 
output fields are coherent. 



It is well known, that the effect of self-phase modula- 
tion due to refractive nonlinearities can be compensated, 
at least in part, by means of an optimum detection pro- 
cedure (for example, by measuring not the phase, but an 
appropriately chosen quadrature amplitude of the probe 
electromagnetic wave) and/or by using non-classical light 
p3| , pT[ . The properties of the input quantum state in the 
methods utilizing non-classical light are thereby chosen 
such that after the interaction the probe wave is in the 
coherent or phase-squeezed state. 



e_(i, L) = e_ 
e+(t,L) = e- 



(68) 
(69) 



The sensitivity of the phase measurement would thus be 
determined by shot-noise only, Acj) = 1/y/Jn). 

In the absence of losses, the sensitivity of the detec- 
tion can even be better than the shot-noise limit, if the 
initial state of the field is appropriately chosen |Q . Mak- 
ing use of a SU(2) Lie-group description, Yurke showed 
that the sensitivity of a phase shift measurement in a 



8 



Mach-Zehnder interferometer can approach the so-called 
Heisenberg limit ~ 1 / (n ) , where (n) is the total num- 
ber of registered quanta ■ 

However, in the presence of losses resulting from the 
resonant coupling the noise compensation by means of 
non-classical light is only partial due to unwanted noises 
added by the medium. Taking into account linear losses 
and assuming, that the entrance field is squeezed in the 
way discussed above, we can rewrite the equation for the 
residual noises in the phase as follows: 

rj{z) = 1 — a.{)Z is the z-dependent transmission coeffi- 
cient. The expression indicates, that for small losses in 
the medium, the noise can be almost completely sup- 
pressed. A maximum signal is achieved however when 
r\ <C 1 and thus the use of non-classical light only leads 
to a marginal reduction of the ac-Stark associated noise. 
This is in contrast to the measurement schemes discussed 
in p3| , pT[ which utilize squeezing to improve sensitivity. 
The change of the expression for the ac-Stark associated 
noise leads to a change of the sensitivity factor / accord- 
ing to 



/ — ► / = 



(l-?7)(ln(77- 1 )+?7-l) 
rj ln 4 (f7 -1 ) 



1/4 



(71) 



certain optimum intensity the fundamental signal-to- 
noise ratio attains a maximum value which represents 
the standard quantum limit of optical magnetometer 
based on phase-shift measurements. This quantum limit 
is determined by the dispersion-absorption ratio of the 
atomic medium and the strength of the intensity-phase 
noise coupling. The unique property of EIT is to pro- 
vide a dispersion-absorption ratio which is independent 
of power-broadening and is given by the lifetime of a 
ground-state coherence. The minimum magnetic level 
shift corresponding to the quantum limit of EIT magne- 
tometers can thus be orders of magnitude smaller than 
that of optical pumping devices. 

We have shown that the best candidate to reach the 
standard quantum limit is a magnetometer in Faraday 
configuration, which has been analyzed in detail. In an 
EIT-Faraday magnetometer the signal reduction due to 
power- and ac-Stark broadenings is compensated by large 
densities of the atomic vapor. The influence of classical 
excess noise is completely eliminated due to symmetry 
and there are much less sources for systematic errors. We 
have also shown that the use of non-classical light and dif- 
ferent detection techniques only marginally improves the 
attainable sensitivity since a maximum signal is associ- 
ated with substantial losses in the atomic medium. 
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We have discussed the influence of ac-Stark shifts on 
the sensitivity of optical magnetometers. We have shown 
that these shifts cause a broadening of the relevant reso- 
nances and give rise to additional noise contributions. In 
absorption-type magnetometers, such as OPMs, the ac- 
Stark associated broadening as well as power-broadening 
lead to a reduction of the signal. We have shown that the 
classical part of these effects can be completely compen- 
sated in an EIT magnetometer in Faraday configuration 
where polarization rotation or, equivalently, the relative 
phase shift of two circular components is measured. 

In a magnetometer based on phase measurements ac- 
Stark shifts lead also to a coupling between intensity 
and phase fluctuations. As a result there are addi- 
tional, ac-Stark associated fluctuations which dominate 
over shot noise beyond a critical laser power. For a 
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FIG. 1. Principle of a drive-probe EIT magnetome- 
ter. Strong drive field in 3-level A system (top) leads to 
transparency of probe field and linear dispersion around 
two- photon resonance 5 = (bottom). Lower plot shows x' 
and x" (real an d imaginary part of probe-field susceptibility) 
in arbitrary units characterizing refractive index and absorp- 
tion. Drive-field Rabi-frequency equals natural width of probe 
transition. 




FIG. 2. A-system in Faraday configuration, 7 r are radiative 
(longitudinal) decay rates, 70 the rate of ground-state coher- 
ence decay (transversal decay); A denotes one-photon- and So 
magnetic-field induced two-photon detuning; E± describe left- 
and right-circular polarized field components. Population ex- 
change (longitudinal decay) between ground-state sub-levels 
is disregarded. Also shown are non-resonant couplings to ex- 
cited states \ci) causing ac-Stark shifts. 
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FIG. 3. Schematic drawing of Faraday measurement. Using 
polarizing beam-splitters the output field is decomposed in 
two orthogonal components E\ and E2 45° rotated relative 
to the x — y system. 
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FIG. 4. Logarithm of minimum detectable Zeeman shift 
in arbitrary units as function of logarithm of laser power 
in units of Po = hvoSnA/X 2 ^- Transmission coefficients are 
77 = 1 - a L = 0.8, 0.1, and 0.01. A0/7 = 10 3 . Also shown 
typical behavior of an optical pumping magnetometer (OPM) . 
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